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ABSTRACT 



A numerical model of a circumstellar debris disk is developed and applied to 
observations of the circumstellar dust orbiting fi Pictoris. The model accounts 
for the rates at which dust is produced by collisions among unseen planetesimals, 
and the rate at which dust grains are destroyed due to collisions. The model also 
accounts for the effects of radiation pressure, which is the dominant perturbation 
on the disk's smaller but abundant dust grains. Solving the resulting system 
of rate equations then provides the dust abundances versus grain size and over 
time. Those solutions also provide the dust grains' coUisional hfetime versus grain 
size, and the debris disk's optical depth and surface brightness versus distance 
from the star. Comparison to observations then yields estimates of the unseen 
planetesimal disk's radius, and the rate at which the disk sheds mass due to 
planetesimal grinding. The model can also be used to measure or else constrain 
the dust grain's physical and optical properties, such as the dust grains' strength, 
their light scattering asymmetry parameter, and the grains' efficiency of light 
scattering Qg. 

The model is then applied to optical observations of the edge-on dust disk 
orbiting ^ Pictoris, and good agreement is achieved when the unseen planetesimal 
disk is broad, with 75 < r < 150 AU. If it is assumed that the dust grains are 
bright like Saturn's icy rings {Qs — 0.7), then the cross section of dust in the 
disk is Afi — 1 X 10^° km^ and its mass is ~ 11 lunar masses. In this case 
the planetesimal disk's dust production rate is quite heavy, ~ 9 M0/Myr, 
implying that there is or was a substantial amount of planetesimal mass there, 
at least 110 earth-masses. But if the dust grains are darker than assumed, then 
the planetesimal disk's mass-loss rate and its total mass are heavier. In fact, 
the apparent dearth of any major planets in this region, plus the planetesimal 
disk's heavy mass-loss rate, suggests that the 75 < r < 150 AU zone at ^ Pic 
might be a region of planetesimal destruction, rather than a site of ongoing planet 
formation. 



1. Introduction 

A debris disk is a dusty circumstellar disk that is often found in orbit about a young star. 
It is also suspected that these dust disks might be sites of ongoing planet formation. This 
thinking is motivated by the dust grains' lifetime due to collisions, which is often much shorter 
than the age of the host star. Evidently, a circumstellar debris disk is also being supphed 
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with fresh dust, and colhsions among unseen planetesimals provides a plausible explanation 
for this dust production. And because planetesimals are also the seeds of planets, it is 
conceivable that a debris disk might be forming planets as well. 

But keep in mind that if the collisional grinding in the planetesimal disk is too vigorous, 
then it is possible that the planetesimals might instead grind away before they have a chance 
to assemble into planets. Indeed, models of the early evolution of the outer Solar System 
propose that the early Kuiper Belt, which is a swar m of comets o rbiting beyond Neptune, 
was initially composed of M^b ~ 30 Earth massed (jKenyon 2002 ). which would be enough 
to allow for the formation of two Neptune-class planets. But that model also shows that 
runaway accretion during the next t ~ 500 Myrs only managed to produce a handful of Pluto- 
sized bodies, while the bulk of the planetesimal mass there remained locked in the smaller 
planetesimals. Meanwhile, collisions among the smaller bodies steadily ground much of that 
mass to dust that is then blown out of the system by radiation pressure, which also stalls 
further g rowth. The mass-loss rate implied by this collisional grinding is ~ I/Mkb ~ 10^'^ 
gm/sec (jStern fc Colwelll 119971 : lKenyonll2002l ). Another outer Solar System scenario is the 
Nice model, which requires its primordial Kuiper Belt to persist until t ~ 700 Myrs since 
formation, which is wh en the giant planets suddenly adjust thei r orbits and trigger the Late 
Heavy Bombardment (iGomes et al.l l2005l : iLevison et al.l |2008| ) . Although the dynamical 
history of the Nice model is rather different from other models, it should be noted that this 
and most other models of the outer Solar System's early evolution predict that the primordial 
Kuiper Belt had a quiescent period lasting several hundreds of millions years, during which 
the Belt would have lost tens of Earth- m asses of material due to collisional g rinding and 
blowout of dust due to radiation pressure Jstern fc Colwelll 1 19971 : iKenvonI [2002f l . 



An interesting question is whether the observed circumstellar debris disks have mass-loss 
rates comparable to that predicted for our early Solar System. The answer to that question 
will then provide some guidance as to whether these disks should instead be thought of 
potential sites for planet formation, or else regions of planetesimal destruction. To address 
this issue, the following develops a model that follows the time-evolution of a circumstellar 
debris disk. Sections I2.1H2.4I derives in some detail the model's physics, but readers not 
interested in those details can skip those sections. Sections I2.4.3f[^^^^ examines how the 
debris disk's structure and appearance depends upon the system parameters, and those 
results can also be quickly gleaned by inspecting the figures there. Section [3] then applies 
the model to observations of the debris disk orbiting /3 Pictoris, with Section H] providing a 
summary of the findings. 



^The lKenvoij (j2002i ) model recommends 10 Earth masses in a 6 AU-wide annulus centered at r = 35 AU, 
but this should be multiplied by ^ 3 to account for the Kuiper Belt's full radial width. 
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The debris disk model 



Radiation pressure is th e dominate perturbation on small dust grains orbiting in a 
circumstellar debris disk (c./.. IStrubbe fc Cliiang|2006l ). and the review in Section [211] shows 
how the resulting dust orbits are simple functions of grain size. Dust collision rates are 
derived in Sections l2.2H2.3t and Section 12.41 shows how to use those rates to calculate the 
time evolution of the debris disk's dust abundance. The remaining subsections then illustrate 
how the simulated debris disk varies with the model parameters. 



2.1. dust orbit elements 

An orbiting dust grain is characterized by its size parameter /3, which is the ratio of 
stellar radiation pressure to gravity. For a spherical grain of radius R, the size parameter is 



_ SL^Qrp _ (LA P \ R 
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where and are the star's luminosity and mass, G is the gravitational constant, c 
the speed of light, p is the grain's volume density, and Lq and Mq refer to solar values 



Burns et al 
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The grain's radiation pressure efficiency is Qrp = Qa + — g)Qs 



19791 ). where Qa is the efficiency of the grain's absorption of starlight and Qs 



is the efficiency of light scattering by the grain. The scattering asymmetry parameter g is 

g= $(0)cOS0dl] (2) 



where the integral runs over all solid angles, and the phase function $((/>) gives the proportion 
of light that is scattered through scattering angle 0, which is the angle between the directions 
followed by the incident and scattered light. The phase function is normalized such that 
J = 1, so [(^l < 1. Forward scattering dust grains have values of g > 0, backscattering 

grains have g < 0, and isotropic light scattering has g = 0. The model developed here will 
be applied to the /3 Pictoris debris disk, whose dust grains are substantially larger than the 
wavelength of the incident light. This is the geometric optics limit, and energy conservation 
in this limit requires Qa + Qs = ^, so 

Qrp = 1 - gQs- (3) 

Evidently, Qrp — 1 except when \gQs\ is not small, which can occur if the dust grains are 
efficient forward or back scatters. Equation ([1]) shows that radiation pressure is significant 
for grains of radii R ~ fim when orbiting a solar-type star, and is unimportant for 



- 5 - 



grains with R ^ 1 /im. Since radiation pressure has the same inverse-square force law 
as stellar gravity, the effects of radiation pressure are easily accounted for by substituting 
GMj, —7- (1 — /3)GM^ into all equations for the dust grains' motion. 

Dust is manufactured when planetesimals collide. Consider a dust grain that is gen- 
erated by colliding planetesimals that reside in nearly circular orbits of radius r. These 
planetesimals have a specific energy Ep = kinetic -|- potential energy that is 

E„ = -v^ = 4 

where v = \fGM^Jr is the planetesimal's velocity, assuming its orbit is circular or nearly 
so. However, a dust grain that forms via a collision at radius r will have a specific energy 
Eld = — GMi,{\ — l3)/r = —GMi,{l — (3) /2a, where a is the dust grain's semimajor axis. 
This assumes that the dust grain has the same velocity v as its parent planetesimal at the 
moment of its creation, or equivalently, the velocity at which the grain is ejected from the 
planetesimal is small compared to its orbital speed. But this is reasonable since the debris 
disks considered here are thin, which implies that dust generation occurs with low ejection 
speeds. In this case, ^v"^ = GM^/2r in the above, so the dust grain's semimajor axis is 



flThebault et al.ll2003[ ) 



= ^zrj^r. (5) 

Bound orbits have a > 0, so only dust having (3 < (3max where f3max = \ will populate the 
resulting debris disk, while dust having {3 > f3max are not bound to the star and quickly leave 
the system. According to Equation ([1]), the smallest bound grains have radii 

Next, note that a planetesimal's specific angular momentum Lp = rv also equals the dust 
grain's angular momentum Ld, since both have the same position and velocity at the moment 
of the grain's creation. Since a dust grain's Keplerian orbit has = a/GM*(1 — /3)a(l — e^) 
while the planetesimal's Lp = \fGM^ = L^,, this provides the grain's eccentricity 

em = ^ (7) 



(IThebault et al.l 120031 ). The dust grain's periapse g^erj and apoapse Qapo distances are then 



Qperi = a{l- e) = r (8a) 
and Qapo = a{l + e) = r/(l - 2/3) (8b) 
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where r is the orbital radii of the planetesimals that gave birth to the dust. Also note that 
periapse is precisely at the planetesimal's orbit, which tells us that the grain's longitude of 
periapse Cj is also the longitude where the grain formed. 



2.2. streamlines 

Assessing the collisions that a dust grain experiences over its lifetime can be a labo- 
rious calculation, one that usually requires evaluating integrals over the dust grains' size 
distribution and the debris disk's three dimensional volume. However those integrals can 
be replaced by a s imple summation when the model disk is quantized; see for example 



Krivov et al.l ( l2006l ). In the following, the planetesimal disk, which is the source of this dust, 
is represented by Nr concentric rings that have radii r uniformly distributed over the inter- 
val Tin < r < Tout- Each planetesimal ring also produces dust at Ni discreet sites that are 
distributed uniformly in longitude about the ring. The dust grains' radii R = {[imax/ P)Rmin 
are also quantized so that they uniformly sample the interval Rmin < R < Rmax, with A^^ 
size-bins in the interval < (3 < f3max- As a result, the planetesimal disk produces dust 
whose grains inhabit Ng = Nj-NiNp distinct streamlines, or orbits, whose shapes are given 
by Equations ([5]) and ([7]), and whose orientations are uniformly distributed in longitude. 

The rate of dust production p{R, r) by the planetesimal disk is also a power law, with 
p{R, r) oc -R"'' being the rate at which the planetesimals inject dust of radii R into a single 
streamline, with q being the power law for the grains' differential size distribution. The model 
will also allow for the dust production rate to vary with the source planetesimals' radial 
distance r as p{R, r) oc {r /routY- A planetesimal disk having c < can then be thought of as 
experiencing 'inside-out' erosion, since the inner portion of that disk experiences greater dust 
production than its outer parts, while a disk having c > suffers 'outside-in' erosion. Also 
note that R oc so the dust production rate can be written r) = po{/3 / PmaxY {r / VoutY 
where po is a constant. 

Since there are Ns = N^NiNf^ streamlines in the model debris disk, the index i will 
be used to identify any one such dusty streamline, where 1 < i < Ng. The streamlines' 
longitudes of periapse Ui are uniformly distributed such that the resulting debris disk has 
A'^-fold symmetry. The quantity ni{t) will be the number of dust grains inhabiting streamline 
i at time t; these grains all have the same size /3j and orbit elements ai,ei,Ui, and they all 
formed at the same longitude Ui in the same planetesimal ring that has radius r,. The 
streamline's abundance ni{t) will then evolve over time t due to the production of dust by 
the planetesimal ring r,, but that quantity will also vary due to collisions with the dust 
grains that inhabit the disk's other streamlines. The rate at which these dust grains collide 
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Fig. 1. — A dust grain of radius Ri crosses streamline j that contains dust of radii Rj. 
Grain i sweeps a length Ai = 2(i?j + Rj)/ sin^jj^ along streamline j where sin^jj^ is the 
angle between the grains' velocity vectors Vj and v^. 

and destroy each other is derived below. 



2.3. collision rates 

Figure [T] shows a dust grain in streamline i as it crosses streamline j. The dust in 
these streamlines have radii Ri and Rj and velocities Vj and Vj. Suppose for now that the 
debris disk is flat, and that all dust grains have zero inclinations. In this case, grain i would 
sweep across a length A£ = 2{Ri + Rj)/ siiKpijk in streamline j, where (pijk is the angle 
between vectors Vj and Vj, so sin (pij^ = |vj x Vj|/t>jt>j. This quantity depends not just on the 
streamlines' orbit elements a^, Cj and aj, Cj, but also on the streamlines' relative longitude of 
periapse Uk = — ^i, which is why the k subscript is also introduced. 
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The linear density of dust grains in streamline j is Xj, so particle i will collide with 
XjAi particles as it traverses streamline j. The flux of dust in this streamline is XjVj, 
which is the rate at which grains cross a point in that streamline, so XjVj = nj/Tj where 

rij is the total number of grains in that streamline and Tj = 27r^a^/(l — f3j)GM^ is that 
streamline's orbit period. So the number of collisions that grain i suffers as it crosses this 
site in streamline j is Ariijk = XjAi = 2(i?j + Rj)njfz/ sm. (f)ijkVjTj, where the factor is 
now introduced to account for the dust disk's vertical thickness. Note that fz would be unity 
if streamlines in the dust-disk where coplanar. However a real disk has a vertical half-width 
h = Irijk, where / is the grains' characteristic inclination and rjj^ is the radial distance 
from the star where streamlines i and j cross, with the k subscript again indicating that 
that this quantity also depends on the streamlines' relative longitude of periapse Uk- Thus 
fz = 2fc{Ri + Rj)/2h = fc{Ri + Rj)/ Ivijk is the probability that two grains in streamlines i 
and j are close enough in the vertical direction to come in contact, and the factor = 7r/4 
is a geometric correction that accounts for the grain's circular cross-section. And since there 
are Ui grains in streamline they will suffer UiAnijk collisions with streamline j after one 
orbit period T,, so the grains in streamline i collide with the grains in streamline j at site Vijk 
at the rate Ui Ariijk /Ti. Thus the total rate at which grains in streamline i suffer collisions 
with all of the disk's other streamlines is 

K. = E ^ = E E E (9) 

j=l k=l ri 

where the coefficient aijk = n{Ri+Rj)^Toutfijk/'^I ^^^(pijkrijkVjTiTj, andTout = 2ti ^/rl^JGMl 
is the orbit period of the outermost planetesimal ring at r = Tout- Note that an additional 
factor f^j^ was also introduced into the above; it takes values of f^^f^ = 1 if the collision with 
the impacting grain Rj is energetic enough to disrupt the target grain and with f^^^ = 
if the collision does not disrupt grain R^. The threshold for coUisional disruption is given 
later in Section 12.3.21 

The leftmost sum in the above runs over streamlines that are composed of dust of radii 
Rj that have size parameters /3j = (3max{Rmin/ Rj) produced by the planetesimal ring that 
has radius r^-; their orbit elements aj and Cj are given by Equations (E} and (j?!). These dust 
grains have Njs possible sizes, and there are Nr such planetesimal rings, so 1 < j < NpN^. 
The middle sum then accounts for the Ni sites in each planetesimal ring that produce dust 
having relative longitudes cD^ = a)j— a)j, so 1 < k < Ni. Also keep in mind that pairs of 
orbits cross at two sites Vijk = ri and rijk = r2 whose contributions are accounted for by the 
rightmost sum. 

But first a comment on differently-sized dust grains that are produced at the same site 
within the same planetesimal ring. Such grains have the same longitude of periapse, Ui = Uj, 
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and thus will re-encounter each other at periapse with the same velocity Vj = Vj. This also 
makes siiKpij^ = 0, which might appear problematic because the collision probability Oijk 
would seem to be singular there. But keep in mind that grains i and j have the same velocity 
at periapse, namely, the planetesimal ring's velocity. Consequently, these grains have zero 
relative velocity at this particular site, so there is no chance for collisional fragmentation, 
and aijk is set to zero in this instance. 

A dust grain's velocity is Vj = {2'Krout/Tout)\/ 0- - l^j){'^rout/rijk - rout/aj) where Tijk is 
the distance from the star where streamhnes i and j cross. Since Ti/Tout = {cii/^outY^'^ / — A, 
the aijk in the above then becomes 

^ out \ I ^out \ I ^min \ /i n\ 

4Jsin0ijfcy 2aj/rijfc - 1 V PiPj ) Krijkaj J V / V f^out J 

This quantity is the collision probability density, in the sense that a^jkUj is the probability 
per time Tout thai a dust grain in streamline i suffers a collision with the dust in streamline 
j at site Tijk. 



2.3.1. orbit crossing sites 

The Vijk in Equation (ITOl) is one of two sites where orbit i crosses orbit j, and it is 
sensitive to the orbits' relative longitude of periapse Uk = ojj — Coi. To solve for rj^fc, rotate 
the coordinate system so that the x axis points to orbit i's periapse. Requiring the two 
orbital ellipses to intersect yields 

'"''^ 1 + Cj cos 6 1 + cos(^^ — Qk) ' 

where Pi = aj(l — Cj) is the orbit's semilatus rectum and 6 is measured from the x direction. 
Using cos(^ — Uk) = cos 6 cos Uk + sin 6 sin Uk allows Equation (ITT]) to be written as A cos 6 — 
Bsm9 + C = 0, where A = e-i —pej cosuk, B = pej sinuk, C = 1 —p, and p = Pi/pj. Solving 
for sin 6*^ = 1 - cos 6*^ then yields cos + 2Dcos9 + E = where D = AC / {A^ + B'^) and 
E = {C^ - B^)/{A^ + B^), which has solution 



cos^ = -D ± VD"^ - E (12a) 
and sin^ = {AcosO + C)/B. (12b) 

Solve Equations (fT2l) for the two 9's and insert those into Vijk = Pi/iX + Cjcos^^) to get 
the radial distances for where these orbits cross. Then rotate the coordinate system back, 
6 ^ 6 + Ui, to get the longitudes where these orbits cross. Note that Equations f[T2|) are real 
when > E, which is the requirement for orbits i and j to cross. 



To evaluate the sin^jjfc in Equation ffTOj) . the dust grain's radial velocity Vi^r and tan- 
gential velocity Vi^ need to be evaluated at the collision site. Those velocity components are 



and Vifi 




= [1 + ei cos(6' - cui)] 



sin(^^ — 



(13b) 



(13a) 



where fij = Iri jTi is the grains' mean motion. These velocity components are then used to 
obtain sin^jj^ = |vj x Vj|/t>jt>j = \vi^r'Vjfi — Vi^0Vj^r\/viVj where fj is the grain's total velocity. 

Note that the above ignored the grain's small vertical velocity Vi^z-, which is of order Ivi. 
This is appropriate since J ^ 1. However the following Section will still need an estimate of 
the dust grains' relative velocity in the vertical direction when determining whether collisions 
are destructive. And since the model being developed here does not faithfully follow the dust 
grains' vertical motions, it is assumed here that any two pairs of dust grains have a relative 
vertical speed that is equal to their rms value, \vi^z — Vj,z\ = \/ (Ivi)^ + {IvjY. 



Now derive the threshold for a dust grain's coUisional disruption. The target dust grain 
has mass mi and velocity Vj, and it is struck by an impacting dust grain of mass mj and 
velocity Vj. The two dust grains' relative velocity just prior to impact is = — Vj, and 
the specific work that the impactor must do on the target grain in order to disrupt it is 
Q* . For simplicity, assume that any debris produced by the disruption of the target grain 
has zero dispersion velocity, which then provides the minimum energy that is needed for 
disruption. In the target dust grain's rest frame, the system's kinetic energy is \mjvl just 
prior to impact. It is shown below that most collisions will shatter both the target and the 
impactor, so the post-impact debris will have a total energy E = ^mjv"^ — Q*{mi + rrij) 
where Q*mi is the work that must do in order to shatter grain rrii. But this energy can also 
be written a.s E = ^rriiv'^ + ^rrijv'^^ where v'^ is the post-impact speed of the debris from 
grain mj, and Vj is the post-impact speed of mass mj in this reference frame. Momentum 
conservation also requires that rrijVr = miv[ + rrijVj. Solving the energy and momentum 
equations simultaneously then provides the post-impact speed of debris m^. 



2.3.2. disruption threshold 




(14) 
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Keep in mind that this is the speed of debris in the reference frame where mass rrii was 
stationary just before impact. Consequently, the speed of that debris in the inertial reference 
frame is 

\i+vlVr (15) 

where is the unit vector that points in the direction of v^. 

The colhsion is energetic enough to disrupt grain rrii when Equation f|T4|) is real, which 
requires > 2(mj + nijYQ* /miirij. Assuming all grains are spheres of similar density, this 
requirement becomes 



vl > 



mi 



(16) 



When the above is satisfied, the collision is fast enough to destroy the target grain which 
is accounted for in Equation f lTOj) by setting ff-^ = 1, while f^,j^ = when the above is not 



satisfied. In that case grain rrii will recoil from the collision and enter a new orbit about the 
star. However, nearly all collisions result in the disruption of grain rrii, so it is safe to ignore 
these rare non-disruptive collision events. 

Unfortunately, the dust grains' collisiona l specific energy Q* i s rather uncertain, and 
likely depends on the target grain's size Ri (IHolsapple et al.ll2002l ). Another d ifficulty is 



that experimental measurements of Q* only extend down to cm-sized target s (IHolsapple 
19941). Howev er models of observ e d deb ris disks suggest ~ 3 x 10^ ergs/gm (IWyatt et al. 
2007bl ). while IStrubbe &: ChiangI (120061 ) recommend a nominal value of Q* ~ 10^ ergs/gm if 
the dust grai ns are rocky. But if cir c umstellar dust grain s are i cy, then values ~ 100 smaller 



are possible ( IHolsapple et al.ll2002l ). IStrubbe fc ChiangI ( 120061 ) also note that sandblasting 
machines accelerate destructive particles to f ~ 100 m/sec, which would suggest an upper 



limit of Q* 



;2 < 



Iqs erg/gm. However, all collisions in this model are disruptive when 



Q* < 10^ ergs/gm, so the following will only consider the interval 10^ < Q* < 10^ ergs/gm. 
And in the example model of Section I2.4.3[ 99.7% of all collisions are destructive when 
Q* = 10'' ergs/gm while 89% of collisions are destructive when Q* = 10^ ergs/gm. 



2.4. the debris disk's time evolution 

The debris disk's dust abundance evolves due to dust production by the planetesimal 
rings minus losses due to collisions among dust grains. To quantify this, let pi = p{/3i,ri) be 
the rate that one site in planetesimal ring rj produces dust of size /3j that gets injected into 
streamline i. The abundance Ui of dust in streamline i then evolves according to the rate 
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equation 




cmt 



Nl T2 

j=l k=\ T\ 



(17) 



where the left term accounts for dust production and the right term accounts for colhsional 
destruction, Equation ([9]). Index i refers to the dust that reside in the target streamhne i, 
while j refers to dust grains in the impacting streamline j that have a relative longitude of 
periapse cj^ = tu-,- — a)j. Although the model consists of Ng = NrNiNjs streamlines, there are 
only Nri3 = NrNjs equations in Equation f[T^ that are distinct, since the A^; streamlines that 
are generated in the same planetesimal ring with the same dust size /3j have an identical 
evolution due to the system's azimuthal symmetry. 

The total number of dust grains of size /3j that are generated by planetesimal ring is 
Ni = NiUi, where Ni is the number of dust-producing sites in a planetesimal ring. Similarly, 
the total rate at which ring rj produces dust of size is Pj = NiPi. Multiplying Equation 
f|T7|) by Nl and setting 



Here, aij is the mean probability per time Tout that a grain of size /3j and orbit elements Oj, Cj 
collides with a grain of size and orbit (3j , aj , ej . Note that this is a mean probability since 
Equation (fT8|) averages probabilities over all possible longitudes of periapse = uj — cD,. 

This model is also going to assume that, when dust grains collide, the resulting dust 
fragments are so small that they are unbound and driven away by radiation pressure. This 
assumption simplifies the problem enormously since, if it where not true, then the number of 
streamlines in the debris disk would grow exponentially as collisions beget additional dusty 
streamlines that collide with other dust and generate even more streamlines. The validity of 
this assumption is confirmed later in Section [3.3.11 




(18) 



k=l ri 



then provides the rate at which Ni evolves over time. 




(19) 



2.4.1. two -component model 



Now consider a rather simple system, that of a single planetesimal ring that manufac- 
tures dust having only two sizes, small (S) grains having radii near the blow-out radius Rmin, 
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and large (L) grains of radii Rl ^ Rmin- In the resulting debris disk, nearly all collisions 
are with the small grains, due to their much greater abundance in the disk, Ns ^ N^, and 
their much greater production rates Ps ^ Pl- Because small grains are only colliding with 
other small grains, their abundance varies as 



dt " T 



out 



according to Equation ( |T9l) . where ass is the probability density for collisions among small 
grains. The solution is 



Nsit) = NlHanhit/ts), (21) 



where 



N'^= and ts= \h^^=^. (22) 

^ ass \ assPs Ps ^ ' 

Here, N'^g is the total number of small dust grains that result when the disk settles into 
equilibrium, which occurs at time t ^ ts when ta.nh.(t /ts) 1 and dNs/dt — )• 0. Note that 
ts is also the small dust grains' collisional lifetime. Evidently the small grains' equilibrium 
abundance as well as their optical depth will vary as N^'^ oc where Ps is their production 

— 1/2 

rate, while the timescale for the disk to settle into collisional equilibrium varies as ^5 oc 

The large grains' abundance also evolves according to Equation (11 9p . which becomes 
dNi/dt = Pl — O-lsNlNs /Tout when collisions with small grains dominate. And when the 
disk has settled into collisional equilibrium, dNi/dt = 0, so the equilibrium abundance of 
large grains is 

Nl" = (23) 

where ais is the large grain's collisional probability density. 

Although this two-component treatment might seem too simple for quantitative work, 
it is illustrative, because a real disk's optical depth and collision rates are often dominated 
by the disk's smallest grains that have size parameters close to Pmax- In that case, the disk's 
collisional equilibrium timescale and its optical depth would all vary with the square root 
of the planetesimal's dust production rate. Equations (122] - [23|) are also useful because they 
provide a convenient test of the more general debris disk model that is developed below 
in Section I2.4.2( using that more general model to simulate a two-component debris yields 
results that agree with Equations ( ETH2^ to within 0.04%. 
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2.4-2. scale-invariant evolution 

Now consider a more realistic scenario where the system's planetesimal rings produce 
dust whose rates are power-laws in the dust size parameter Pi and planetesimal ring radius 
Tj such that Pi = Poi^Pi/ [3maxYi.fi/routY- It will also be convenient to first convert Equation 
ffT^ into a dimensionless system of equations via the substitutions 

t = Toe and N,{t) = Nj^\ Ntif) (24) 

\ Pmax J \ Tout J 

where t* = t/To is dimensionless time coordinate and N*(t*) = Ni(t)Po/PiNQ is the scaled 
abundance of grains of size Pi produced by ring r^. The collision probability density will 
also be written as aij = {a*^/ I){Pj/ Prnax)~'^ijj/rout)~''{Rminl'f'outY where the scaled collision 
probability is 



sin^jj-fcU 2aj/rijk — I \ PiPj / \ Pmax J K^ijkO'j J \ C'i J \f^out J 

In the above, the constant timescale To and population scale-factor A'o are 



Pq \ R, 



To = jll^fll^] (26a) 



and No = P,To=^/IP^t(^)- (26b) 
Inserting this into Equation f[T^ then yields a system of equations 



dm 



NtY^a^N* (27) 



dt* 

that is not just dimensionless, but also scale invariant. This is very handy, since Equation 
( |27|) need only be solved once in order to apply the result to a variety of systems that might 
have a range of planetesimal radii rout orbiting stars of varied masses and luminosities L^. 
Note, however, that the rescaled collision probability a*- does depend of the dust power laws 
q and c, the dust grains' specific energy for collisional disruption Q*, and the planetesimal 
disk's radial width r in /rout (Equation [25]), so any changes to those parameters does require 
solving Equation (127|) again. 

Equation fl27|) is a coupled set of nonlinear differential equations whose initial conditions 
are N*{0) = 0. These equations are easily solved numerically for each streamline's relative 
abundance N*{t*) using a Runga-Kutta algorithm. That solution to Equation f l27|) then 
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provides the dust grains' relative abundances Ni{t)/NQ = {Pi/ Pq)N*{V'). A sample calcula- 
tion is shown in Figure |2l which plots the relative abundances N^/Nq versus dimensionless 
time t* for dust that is generated by a narrow planetesimal ring (A^^ = 1 and rj„ = Tout) 
whose dust production has the q = 3.5 size distribution that expected for debris that is in 
collisional equilibrium ( iDohnanyil Il969l ) . In this and all calculations that follow, the simu- 
lated debris disk is composed of grains having = 200 distinct sizes, with the dust size 
parameters distributed over the interval 0.0652 < f3i < 0.497 such that the dust grains' radii 
Ri = {Pmaxl Pi)Rmin are Uniformly sampled over the interval Rmin < -R < Rmax- Dust of each 
size Pi are launched from Ni = 100 evenly-spaced longitudes about the planetesimal disk, 
so the resulting debris disk is comprised of A^^ = A^A^jA^ = 2 x lO^A^^ distinct streamlines. 
The planetesimal disk that is the source of this dust is usually composed of A^r = 3 to 5 
planetesimal rings, except A^,. = 1 is used when the dust source is a narrow planetesimal ring, 
as is the case for Figure [2J The /3 distribution used here results in streamlines that range as 
far as r = IdOvout from the central star, where Tout is the planetesimal disk's outer radius. 
Execution times on a desktop PC range from 4 minutes for simulations having A^ = 1 to 2 
hours when A^ = 5. 
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Fig. 2. — Equation fl27|) is solved for the scaled dust abundances N*(t*) assuming that the 
dust is produced by a narrow planetesimal ring {Nr = 1) with a. q = 3.5 size distribution. 
The scaled dust abundances are then converted into absolute dust abundances Ni(t*)/NQ 
and plotted versus dimensionless time t* (left and lower axes). This model also assumed 
that the dust are weak, Q* = 10^ ergs/gm, so all collisions disrupt the dust. Abundances for 
nine selected streamlines having the indicated size parameter /3 are shown above; not shown 
are the model's 191 other curves that behave similarly. The right and upper axes are scaled 
for an Tout = 50 AU planetesimal ring in orbit about a solar-type star that produces dust at 
the rate = 10^^ gm/sec with p = 1 gm/cm^, Q^p = 0.96, and / = 0.1 radians, for which 
Nq = 2.8 X 10^^ grains and Tq = 5.6 x 10"^ yrs. The dashed grey curve is Equation f ETj) with 
an equilibrium timescale T^q = 7.1 Myrs. 

Although the calculation seen in Figure H] covers a fairly broad range of (3 parameters, 
the resulting debris disk is nonetheless populated by dust having a fairly narrow range of 
sizes, since it is only the smallish dust grains that are lofted well outside of the birth ring 
by radiation pressure. In Figure [21 the central star is solar and the dust grains have a 
density p = I gm/cm^ and radiation pressure efficiency Qrp = 0.96, so the radius of the 
smallest bound dust grain is Rmin = 1-1 pm, and the range of dust radii that populate 
the simulated debris disk is only 1.1 /im< i? < 8.5 pm. Of course, a real disk will also 
manufacture larger dust, but those larger grains would be confined to the disk's innermost 
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region r^n < r < l.lhrout, so they do not alter the disk's large-scale structure at r ^ Tout- 
Also, these larger grains contribute little to the disk's optical depth due to their slower 
dust production rates and their very short collisional lifetime, so their absence from these 
simulations is justified. 

The parameter Pq that appears in the constants A'^o and Tq (Equations [26i) is related 
to the planetesimal disk's total dust mass production rate, which is Md = J^i^i^P^i 
where the sum proceeds over all streamlines that make up the dust disk. Since Pi = 
Po{(3i/ (^maxYiri/routY, this provides Pq = Ma/ Sgcmmin, where rrimin = ^pRlain is the mass 
of the smallest bound dust grain, and the factor Sqc = t^maxY'^iri/routY- Inserting 

this and Equation ([6]) into (1261) then provides the scale factors 



P f rout f Sqc'^ ""^/^ 



Igm/cm^; V50AU/ V107 



(28b) 

and To = 5.6 x lO^Q^/^ ( ] ( — i— ) ( ) ( ) (28c) 

1013 gm/sec^ VO-lrady \LqJ \MqJ ^ ' 

( Tout y/u^cicY^^ 

nsOAu) (,107 J ^'"^^ 

where the factor Sqc = 107, 320, and 530 for the simulations having A^^ = 1, 3, and 
5 planetesimal rings with c = and q = 3.5. The above quantities allow one to easily 
rescale all figures shown here for systems having alternate dust production rate or the 
planetesimal disk's outer radius Tout, etc. The following subsection will also show that the 
time Teq for the resulting dust-disk to settle into an equilibrium where dust production 
balances mass-loss due to collisions is T^g ~ 130To. 



2.4-3. example: ar = 50 AU birth ring 

To illustrate, consider a narrow (A^ = 1) planetesimal ring of radius rout = 50 AU in 
orbit about a solar-type star. This ring will have a dust mass production rate of M^i = 10^3 
gm/sec, with dust having a density p = 1 gm/cm3 and inclinations / = 0.1 radians. The 
grains are dark, with Qs = 0.1, and are asymmetric light scatters having g = 0.4, so the 
radiation pressure efficiency is Qrp = 0.96, A'o = 2.8 x 103^ and Tq = 5.6 x 10^ yrs. Multiplying 
the lower and left axes in Figure |2] by these constants then provides the dust abundances 
Ni{t) versus physical time t, which can also be read off the right and upper axes of Figure 
[21 That figure shows that the smallest dust grains having /3 ~ Pmax need the most time 
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Fig. 3. — The black N{R) curve is the differential size distribution for the dust in the debris- 
disk model of Figure [2] after the disk has achieved equilibrium at times t ^ T^q assuming the 
dust grains are weak with Q* < 10^ ergs/gm. Grey curves are for models having stronger 
dust grains, Q* = 10'' and 10® ergs/gm. The dashes indicates an size distribution, and 
the dotted curve is proportional to the planetesimal ring's differential dust production rate 
P{R) oc R-^-^. 



to settle into an equilibrium where dust production balances destruction due to collisions. 
The grey dashed curve in Figure |2] also shows that the abundance of the smallest dust does 
resembles Equation fl^T]) when the equilibrium timescale T^q = 7.1 x 10^ yrs is chosen as the 
moment when the smallest dust have reached two-thirds their equilibrium abundance. This 
demonstrates that the two-component dust size model of Section 12.4.11 is in fact relevant to 
this kind of debris disk. And because Tq and T^q have the same dependance on the model 
parameters, comparison to Equation ( I28cp shows that T^q ~ ISOTq. 

After a time t ^ Teq, the debris disk will have settled into equilibrium. The solid 
black curve in Figure [3] shows the disk's equilibrium grain size distribution N{R), which is 
extracted from the rightmost part of Figure [2] and plotted versus grain radii R. That model 
assumes that the dust grains are weak, Q* < 10^ ergs/gm, which means that all collisions are 
destructive. Note that the debris disk's dust size distribution N{R) is significantly steeper 
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than the q = 3.5 size distribution that governs the ring's dust production (the dotted curve 
in Figure [3]), with A^(-R) oc except near R ~ Rmin, where it has an even steeper 
deperi dance. A peak at ~ Rmin is also seen in the debris disk modej^ of Krivov et aL 



( 120061 ) ■ Figure [3] illustrates the main consequence of dust-dust collisions, which tends to 
destroy the disk's larger dust grains at a faster pace. However, when the dust grains are 
stronger, with 10'' < Q* < 10^ ergs/gm, then the larger grains are more resistant to collisional 
destruction, and they become more abundance (Figure [3l grey curves). 

When the system is in equilibrium, the rate at which the planetesimal ring injects 
dust of size (3i into the debris disk. Pi, balances the rate at which collisions remove dust 
from the disk TZi, so TZi = Pi. Thus the collisional lifetime of the grains in streamline i is 
Tc{Ri) = Ni/TZi = Nil Pi, which is plotted versus grain radius R in Figure S] (see the black 
curves that are labeled by their dust production rates M(i). That lifetime is simply the time 
when the dust abundances Ni{t) seen in Figure [2] flatten out. As expected, the smallest 
grains that have size parameters /3 that are just shy of Pmax = \ are very long lived. This 
is due to their orbits having very large apoapses Qapo = fout/i)- — 2/3) ^ Vouti so small dust 
spend most of their time far from the planetesimal ring, in regions where the disk's optical 
depth is low and collisions are rare. Note also that Tc oc (A^o/-Po) oc Tq, so the dust grains' 

. 1/2 

collision timescale obeys the same scaling as Equation fl28cp . in particular with oc 
Figure m also shows that nearly all dust grains in such a debris disk have lifetime s that are 



short compared to the age of the host star, which is typically ~ 10'' to 10^ years (jMeyer et al. 



2003). 

When the grains are weak, Q* < 10® ergs/gm, then the lifetime of the larger R > 
2Rmin varies at Tc{R) oc R~'^-^ (Figure Hj). That figure also shows that increasing the dust 
grains' strength Q* increases the larger grains' longevity, without affecting the small grains' 
collisional lifetime. This is because all collisions are energetic enough to destroy all small 
grains. 



Note though that the equilibrium dust size distribution in the iKrivov et aD |2006l ) model is 'wavy'. 
Waves in a size-distributi on occur when som e process tends to fa vor the rapid removal of the smallest bodies 
(jPurda fc DermottI 119971 ) ■ and Figure 5 of iKrivov et al.l (|2006f ) would suggest that next peak in our size 
distribution might occur at sizes R ^ 1 mm. However these large grains would still have a very short 
collisional lifetime (see Figure HJ, which makes their influence in the dust-disk quite negligible. 
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Fig. 4. — The lower middle black curve is the dust grain's collisional lifetime Tc{R) plotted 
versus dust radius R for the debris-disk model of Figure [2] that has a dust production rate 
Md = 10^^ gm/sec and strength Q* < 10^ ergs/gm, while the other black curves are for 
systems that are otherwise identical but with the indicated dust production rates M^. The 
grey curves are for stronger dust grains having Q* = 10^ and 10* ergs/gm and Md = 10^^ 
gm/sec; compare also to Figure [31 The upper solid grey curve Tq is the timescale over which 
the dust grains' semimajor decays due to PR drag, with Tg the timescale for eccentricity 
damping due to PR drag; see Equations (1361) . 



2.4-4- optical depth 

To determine the model debris disk's normal optical depth x„(r), it is convenient to first 
calculate the total dust cross section A[r) that resides interior to distance r from the star. 
That quantity is obtained by first counting the number of dust grains Anj{r) in streamline 
j that also lie interior r, which is Anj{r) = J Xjdi where the integration runs along the 
stretch of streamline interior to r. This becomes a trivial integral over time t after noting 
Xjdi = XjVjdt where Xj = rij/vjTj is the dust grain's linear density, and rij is the total number 
of grains in streamline j that have velocity Vj and orbit period Tj. Consequently, Anj(r) = 
2njtj{r) /Tj where tj{r) is the time for dust in streamline j to travel from periapse to distance 
r. And since rii = Ni/Ni where Ni is the total number of dust grains in the disk that have 



10 



distonce r 




Fig. 5. — The dust optical depth r„(r) is plotted versus distance r in units of the planetesimal 
disk's outer radius Tout for the model of Figure [2l Shown is r„(r) at dimensionless times 
t* = 0.01, 0.1, 1, 100, and an r'^/^ curve. 

the same sizes and orbits, the total cross section for grains of radii Ri = {PmaxI l^ijRmin that 
reside interior to r is Ai{r) = NiAni{r)TTR'f = 2Ni{ti/Ti){Praaxl PifT^Rl^in- The total dust 
cross section interior to r is then A{r) = Ai{r) where the sum proceeds over all streamlines 
whose dust have distinct sizes /3j produced by planetesimal ring r^. The ratio tj/Tj in the 
above is obtained by solving r = aj(l — Cj cosE) for the dust grains' eccentric anomaly E{r), 
which is then inserted into Kepler's equation to obtain tj/Tj = {E — Ci sin£')/27r. Also note 
that the differential AA = {dA/dr)Ar is the dust cross section that resides in an annulus 
of radius r and width Ar. Since the disk's normal optical depth r„ is simply the surface 
density of dust cross section, that quantity is 



which is easily calculated by differentiating A{r) numerically. 

Figure [5] plots the dust optical depth T„(r) versus distance r from the star at dimen- 
sionless times t* = 0.01, 0.1, 1, and 100 for the model of Figure [21 As the figure shows, it 
is the outer portion of the disk that is populated at later times by dust. This is due to the 
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Fig. 6. — Black curves give the debris disk's equilibrium optical depth T„(r) versus radial 
distance r for dust generated in a narrow debris disk of radius a^. Model parameters are 
identical to those adopted in Figure |2] except that different dust production rates are con- 
sidered, Md = 10^^, 10^^, and 10^^ gm/sec, and weak grains having Q* < 10^ ergs/gm, while 
the one grey curve assumes strong dust having Q* = 10^ ergs/gm. 



outer disk being composed of smaller dust that, according to Figure [21 are the last to arrive 
at a collisional balance. Also note that when the disk has settled into equilibrium at times 
t* ^ 1, th e outer disk at r ^ r,„,,t ha s an optical depth r„(r) oc r~'^/^, which is in agreement 



with what IStrubbe &: Chiangl ( 120061 ) call a type B debris disk 



The black curves in Figure [6] show the equilibrium optical depth Tn{r) for simulated 
debris disks that are identical to the one considered in Figure [21 except that different dust 
production rates are considered. These curves assume the dust grains are weak, with 
Q* < 10^ ergs/gm. As expected, these optical depths vary as (see Equation I28ap . 

The one grey curve there shows the optical depth of a disk composed of strong dust having 
Q* = 10^ ergs/gm. That disk shows a slight overdensity near the planetesimal ring, and is 
due to an excess of larger grains in low-eccentricity orbits (see grey curves in Figures [3] and 

mi- 
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Fig. 7. — The normal optical depth r„(r) for three debris disks. One is generated by a narrow 
planetesimal ring of radius Tout (black curve), another by an intermediate-width planetesimal 
ring whose inner radius r^n = OSvout (grey curve), and the third due to a broad planetesimal 
disk having rj„ = O.dvout (dashed curve). The narrow planetesimal disk is represented by 
one planetesimal ring (A^^ = 1), while the intermediate and broad disks use = 3 and 
Nr = 5 planetesimal rings, respectively. All simulations have c = 0, so the dust production 
rate is independent of distance r in the planetesimal disk. Model parameters are otherwise 
identical to those adopted in Figure [2j Also shown is an r~^/^ curve. 

The sharp peak in T„(r) at r ~ Tout seen in Figure E] is due to the narrow width of 
the planetesimal ring that is the source of this dust. However, Figure [7] shows that peak 
broadens when the planetesimal ring's radial width is increased. 

2.4-5. disk surface brightness 

This subsection will calculate the surface brightness of starlight that a simulated debris 
disk will scatter towards an observer that views the disk at optical wavelengths. In the 
following, flux F is the power per area in the incident or scattered radiation, intensity / is 
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Fig. 8. — The observer's line-of-sight (LOS) passes through streamhne j and intersects 
segment 6i at polar coordinates (r, 6'); that segment also lies at a projected distance x from 
the central star at the origin. The insert zooms in on segment 6i, which has a projected 
length 6x when viewed by the observer. The dust in that segment will have traveled a 
distance 6i = Vj6t in time 6t, which corresponds to radial and tangential displacements 
Vj^rSt and Vj^edt, so the angle ip satisfies sin ip = Vj^/vj and cosy? = Vj^r/vj- Also note that 
the scattering angle (p is the angle between the radial direction f and the direction to the 
observer, so = 7r/2 + ^. 



the power per solid angle in a pencil-beam of radiation, and the surface brightness B is the 
radiation beam's power per area per solid angle. 

Begin by calculating the surface brightness of a small segment of length 6ij in streamline 
j. Summing the contributions from all such streamlines will then provide the debris disk's 
total surface brightness along some line-of-sight. Since these disks are often observed nearly 
edge on, an edge-on viewing geometry is also assumed here, though these results are easily 
generalized for other viewing geometries as well. Segment 6ij in streamline j is composed 
of dust having a total cross section 6a j, so the intensity of starlight 61 j that is refiected 
by that segment is 6Ij = Qs^Fi6aj, where Fi = L^/47rr^ is the fiux of incident starlight, 
is the central star's luminosity, and r the dust grains' distance from the star. The dust 
grains' phase function is $, and it is normalized so that its integral over all solid angles is 
unity, with this quantity having units of steradians"^. Note also that Fi6a is the power in 
the incident starlight while J 6IjdQ (when integrating over all solid angles dQ) is the power 
in the scattered light, so Qs = f 6Ijdfl/Fi6aj is the efficiency of light scatter ing by these 



dust grains. In planetary astronomy this quantity is known as the bond albedo (ILester et al. 



19791). 
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The surface brightness of the small segment is Bj = SIj/6Aj where 6Aj is the projected 
area on the sky occupied by that segment whose total length is 6ij and projected length 
(as seen by the observer) is 6xj; see Figure Ei Since that segment is actually a ribbon of 
material of height 2/r due to the dust grain's inclinations /, the segment's projected area is 
SAj = 2IrSxj, and its surface brightness contribution is Bj = {Qs^Fi/2Ir){6aj/6xj). 

The number of dust grains in segment 6ij is 6nj = XjSij = \jVj6t where Vj is the 
dust grains' velocity there and 6t is the time for the dust to traverse 6ij. But this becomes 
6nj = {nj/Tj)6t where Tj is the dust grains' orbit period, since XjVj = rij/Tj where uj 
is the number of dust grains in streamline j that also have radius Rj = {Pmax/ Pj)Rmin- 
Consequently, the cross section of dust in segment Sdj is 5a j = 7inj{l3max/ /3j)'^Rfnini^'^/'^j)■ 
Th.Sit segment has length 6ij = Vj6t, and Figure [H] shows that its projected length is 6xj = 
Vj \ cos{6 + ip)\6t where 6 is the segment's longitude as measured from the x axis. The angle 
ip obeys simp = Vj^e/vj and cosy? = Vj^r/vj (see Figure [8]) where vj^r and vj^e are the dust 
grain's radial and tangential velocities. Equations (fT3l) . so 6xj = \vj^rCos9 — Vj^0sm9\6t = 
2TTaj \ sin 9 + ej sin Cjj\{6t/Tj)/ y^l — . Inserting these results into Bj and noting that rij = 

N,/Ni then yields B, = B,^l - e]{m,N,/Ni){Rl^y^Jr^a,){Pmax/ Pjf /\ sm^ + e, sinc^,| 
where the constant 

Note that fii = 1 steradian is introduced into the above so that Bq has units of surface 
brightness and that the combination <l>f2i is dimensionless. Summing the contributions from 
all streamlines then yields the disk's total surface brightness B{x) = "^Bj as a function of 
projected distance x = r cos 9 from the central star, so 

Bix) - B^l^l^l^ ^— -J [^J[-^ (31) 

1=1 j=l near ' ■' jix / x / \ > / 

where the innermost sum is over the two segments — one nearer and the other further from 
the observer — that intersect the observer's line-of-sight (see Figure E]), while the middle sum 
proceeds over each streamlines' orientation uj, and the leftmost sum proceeds over the Nrfs 
streamlines having distinct dust sizes Pi that originate in the various planetesimal rings rj. 
Equation fl3T]) also requires the two longitude 9{x) where the observer's line-of-sight intercepts 
streamline j. That is obtained by solving x = r cos 9 for 9{x) where r = pj/[l + ej cos{9—u}j)], 
which yields 9{x) = Uj — ta.n~^{B/A) ± cos~^ {x / yA^+B^) where A = pjcosuj — CjX, 
B = pj sinuj, and pj = aj{l — e^) is the ellipse's semi-latus rectum, with the r in Equation 
(pTj) from r = xj cos 6'. 
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Fig. 9. — Equation ( I3T1) is used to calculate the surface brightness B{x) versus projected 
distance from the star x for three edge-on debris disks that are generated by narrow planetes- 
imal rings that have the indicated dust production rates = 10^^, 10^^, and 10^^ gm/sec. 
Black curves assume weak dust having Q* < 10^ ergs/gm, a light scattering asymmetry 
parameter g = 0.4 and albedo Qs = 0.1, with all system parameters identical to those of 
Figure [2l These curves give the disks' unfiltered surface brightnesses integrated over all op- 
tical wavelengths. A x""^/^ curve is also shown, and the grey curve is for strong dust having 
= 10* ergs/gm. 



The phase function employed here is the Henyey-Greenstein function 



$(0) 



ster 



(32) 



47r(l + g^- 2^ cos 0)3/2 

that is widely used in studies of circumstellar dust. The scattering angle in the above is 
related to the dust grain's longitude 9 via = n/2 + 9; see Figure [HI This phase function 
is controlled by the dust grains' asymmetry parameter g = J $(0)cos0(in (see Equation 
E]), with a positive value resulting in forward light scattering while negative values result in 
backscattering. 



Equations f l3TI - [32|) are used to calculate the surface brightness for three edge-on debris 
disks whose parameters are identical to those of Figure |2] except for differing dust production 



u 
u 



o 



E 
u 

m 

CP 

\_ 

0) 



X 

(/) 

CO 

c 

sz 
gi 

(D 
U 
O 

D 



-27- 



10-" 


g = 0.4 ^\ 


- 


10-^2 






g=o ""^^^X 


- 








10-^" 



















0.1 



1.0 



10.0 



projected distonce from star x (r^ ^) 



Fig. 10. — Surface brightness B{x) is plotted versus projected distance x for three edge-on 
debris disks that are generated by a broad (rj„ = 0.5rout) planetesimal disk having a dust 
production rate = 10^^ gm/sec, a hght scattering efficiency of Qscat = 0.1, and the 
indicated hght scattering asymmetry parameter g = 0, 0.4, and 0.8. Also shown is an 
curve. Note also that the small bump in B{x) at x = O.bVout (which is where the line-of-sight 
runs along the planetesimal disk's inner edge) gets washed-out when the light scattering is 
very asymmetric with g > 0.8. 



rates M^; see Figure M As expected, these surface brightness curves vary as i? oc vM^. 
Also note that B( x ) oc y "^^^, which again is in agreement with the type B debris disk of 
Strubbe fc Chianel J2006h . 



Figure [10] illustrates how the inner part of a debris disk's surface brightness depends on 
the degree of asymmetry in the dust grains' light scattering. For instance, when the dust 
grains are either isotropic light scatters {g = 0), or are only modestly forward scattering 
{g = 0.4), the debris disk's surface brightness profile B{x) is roughly constant in the inner 
regions where x < Vout- However, if the light scattering by the dust is strongly asymmetric, 
with g = 0.8, then B{x) continues to increase inwards of x = Tout- Disk s having this kind 



of knee-bend in their s urface brightness profiles appear to be common (IKrist et al.l 12005 



Golimowski et al.ll2006l ). with forward scattering of starlight being the preferred explanation 
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for the bent surface brightness profile of th e edge-on debris disk at AU Mic and /3 Pic 



(IStrubbe &: Chiang! l2006l : lAhmic et al.l |2009| ) . Evidently, light scattering by circumstellar 
dust is rather asymmetric. Lastly, note that almost identical results are obtained when the 
dust grains are backscattering, with g < 0, with those slight differences being due to the dust 
grains' slightly larger radiation pressure efficiency Qrp = 1 — gQs- And the following will 
assume that the grains are forward scattering, which is appropriate for dust grains that are 
larger than the wavelength of the incident r adiation, and is consis tent with forward scattering 



by the Solar System's interplanetary dust (ILamy &: Perrinlll986l ). 



Figure [TT] also shows how an edge-on disk's surface brightness B{x) evolves over time. 
As expected, B{x) gets shallower at later times as the smaller dust grains steadily populate 
the outer parts of the disk; see also Figures |2] and O Eventually, the smaller grains arrive at 
collisional equilibrium at dimensionless times = t/To ^ 1, and the disk surface brightness 
settles into the expected B{x) oc x~^'^ power law. But note that at earlier times, such as at 
time = 0.01, the disk's surface brightness can be as steep as B{x) oc Consequently, 
disks having a surface brightness steeper than x~^'^ might indicate that the disk is younger 
than time Tq, or that dust production by the planetesimal disk has increased in recent times. 

Note that these surface brightness calculations assume that the dusty disk is so tenuous 
that the grains do not shadow each other; this assumption is confirmed in Section 13.3.21 Also 
keep in mind that Figures [914TT] give the disks' total surface brightness integrated over all 
optical wavelengths. If these systems were instead observed through a narrowband filter hav- 
ing a central frequency u then the star's total luminosity = An^Rl B^du = AirRlaT^ 
(where is the stellar radius, B^, is the Planck function, a the Stefan-Boltzmann constant, 
and is the star's effective temperature) in Equation ( 130!) should instead be replaced by 
its specific luminosity dL^/du = nL^Bi^/crT^, so that Equations (13T14321) would then provide 
the disk's surface brightness per unit frequency. 



2.4-6. a relic debris disk 



It is also worth considering the debris disk's evolution when the planetesimal disk's dust 
production suddenly ceases. Such might occur if there are any planets in the system that 
quickly adjust thei r orbits, perhaps d ue to a rapid migration of a planet through a dense 
planet esimal disk (IGomes et al.M2004l ). or due to planet-planet scattering {i.e., the Nice 



model. Gomes et al. 



2OO5I ). Should that occur, then the planetesimals might find themselves 



in unstable orbits, which then leads to a rapid dynamical erosion of the disk and a cessation 
of dust production. 



(N 
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Fig. 11. — Surface brightness B{x) is plotted versus projected distance x (in units of rout) 
for the debris-disk model that is described in Figures |2] and [5] viewed edge-on with Qg = 0.1. 
Shown is B{x) at dimensionless times t* = t/Tg = 0.01, 0.1, 1, and 100, as well as power 
laws that vary as x~^'^ and 



Halting dust production is simulated here by setting the planetesimal ring's dust produc- 
tion rate -Pj(t) = at all times after tstop, or equivalently changing the 1 — )• in the invariant 
evolutions equation ( 127|) at dimensionless times t* > t*^op where t^^^p = tstop/To. Solving those 
equations numerically shows that the larger grains in the inner part of the debris disk are 
quickly destroyed due to their short coUisional lifetimes (Figure S]). This depletes the inner 
disk's optical depth and drives it towards a shallower T„(r) oc r^^/^ power law that decreases 
over time. And as Figure [12] shows, the cessation of dust production also decreases the inner 
part of an edge-on disk's surface brightness so that B{x) oc at times t* ^ t^top- So the 
detection of a debris disk having a shallow optical depth profile T„(r) oc r~^/^ (when the view 
to the disk is oblique or face-on), or an edge-on disk having a shallow B(x) oc x~^^'^ surface 
brightness profile, would indicate that the object is a relic disk wherein dust-production has 
ceased. That relic disk will be composed mostly of small, marginally-bound grains of radii 
R ~ Rmin- The abundance Ns of those small dust grains will then fade over time according 
to Equation (1201) . whose solution is Ns{t') = Ns{tstop)/[^ + N s{t stop) (isst' /Tout] when the 
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The dust-producing debris disk model of Figure [2] is evolved until dimensionless 
= 100 when dust production ceases (Pj = 0). The debris disk's surface brightness 



B(x) is shown at the indicated times t*, as well as two power laws x '^^^ and x 



dust production rate Ps = 0, where t' = t — tstop is the time since the end of dust production 
and ass is the small grains' collisional probability density. Note that this behav ior over time 
is typical for a system tha t steadily grinds away without any replenishment (cf. . IWyatt et al. 
2007al : iLohne et al.ll2008[ ). Inspection of the upper two curves in Figure [12] also shows that 



the time for the system to transform into a relic disk is about ttrans ~ IOTq ~ O.lTgg i.e., the 
transition occurs relatively quickly in comparison to the system's equilibrium timescale Tj 



eg- 



3. Application to /3 Pictoris 



The colored curves in Figure [13] show the surface brightness of the starlight that is 
scattered by the edge-on debr i s disk orbiting /3 Pictoris. These curves are extracted from 
images that iGolimowski et al.l ( l2006l ) acquired with the Hubble Space Telescope (HST) at 
optical wavelengths. These profiles exhibit the classic signature of a debris disk that is 
generated by a disk of colliding planetesimals that extends out to about Tout — 150 AU, 
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with the surface brightness falhng off steeply as B{x) oc x~^'^ at projected distances x > 220 
AU, and less steep interior to Vout- Note also that the /3 Pic disk is asymmetric, with the 
outer part of the disk's northeast (NE, blue curve) ansa being about 50% brighter than 
its southwest (SW, red curve) ansa. Although that asymmetry is not accounted for by the 
axially symmetric debris-disk model that is used here, possible causes for that asymmetry 
are described in Section [3. 3. 5[ 

/3 Pic is an A5V star that lies a distance A = 19.28 pc awa y, has a lurninosit y = 8.7Lq, 
mass ~ l.SMf n, an effective ternperat ure T = 8200 K (jCrifo et al.l 119971 ). and an age 
of ~ 12 Myrs (IZuckerman et al.ll200ll ). Note that the star's peak emission occurs at 
wavelength A = 0.35 /im, yet the radius of the smallest bound grain is Rmin — 3 /xm 
assuming the grains have a density of p = 1 gm/cm^ and a radiation pressure efficiency 
Qrp = 0.53 (Equation [H]). So this system is in the geometr ic optics limit, an d the light- 
scattering theory employed here is appropriate. Figure 11 of iGolimowski et al.l (120061 ) also 
shows that the disk's vertical distribution has a full width at half maximum (FWHM) that 
varies as FWHM oc 0.14x ~ 2h, so the dust grains' inclinations are approximately I = h/x ^ 
0.07 radians. 



The black curve in Figure [13] shows that the debris disk model can reproduce (3 Pic's 
observed surface brightness when the dust-producing planetesimal disk is quite broad, with 
inner and outer radii = 75 AU and Tout = 150 AU. This system has been evolved 
over the star's lifetime, = 12 Myrs, so the system is still approaching equilibrium since 
Teg = 20 Myrs. The dust grains must also be rather asymmetric light scatters, with g = 0.67 
(assuming forward scattering), in orde r to account for the knee seen in the surface brightness 
profile at X ~ 150 AU, consistent with lAhmic et al.l (120091 ). The model also assumes that the 
dust grains are rather reflective, Qg = 0.7, which is the B ond albedo of Saturn's icy A and 
B rings at the observation wavelength (jPorco et al.ll2005l ). With these assumption in hand, 
fitting the simulated surface brightness to the disk's observed B{x) requires the planetesimal 
disk's dust production rate to be Md = 1.7 x 10^^ gm/sec. The total cross section of dust 
in the simulated debris disk is = 1.9 x 10^*^ km^, and the total mass of d ust is = 1 1 



l unar masses. Note that this mass is comparable to the 8 lunar masses that iHoUand et al. 



(119981 ) inferred from submillimeter observations of this disk, which supports the contention 
that the /3 Pic dust grains are rather reflective. 

The dust production rate inferred here is about — 9 M^/ Myr. For comparison, 
this rate is about 200 times higher than what IStern fc Colwelll (119971 ) report for their model 
of the coUisional erosion of our own Kuiper Belt. Note that /3 Pic's dust production rate 
is quite considerable, for if it has held steady over the system's lifetime, then the unseen 
planetesimals orbiting (3 Pic would have lost a total mass of Mdt^ ~ 110 earth-masses. 
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Fig. 13. — Red and blue curves give the surface briglitness B{x) of the edge-on disk 
orbiting Pictori s , extr acted from the Hubble Space Telescope observations acquired by 
Golimowski et al.l ( 120061 ) in the F606W filter, and plotted versus projected distance x from 
the central star. The black dashed curve is for a model that adopts a g = 2.5 size distri- 
bution for dust produced by a broad disk of planetesimals orbiting at 75 < r < 150 AU 
from the central star; this system was evolved for t = 1.2 x 10^ yrs, which is the age of the 
central star. The simulated dust grains' inclinations are / = 0.07 rad = 4.0°, their density 
is p = 1 gm/cm'^, and they are strong with = 10^ ergs/gm. The dust grains also have a 
light scattering efficiency Qg = 0.7, a light-scattering asymmetry parameter g = 0.67, and a 
radiation pressure efficiency Qrp = 1 — gQs = 0.53. The planetesimal disk's dust production 
rate is Md = 1-7 x 10^^ gm/sec, so Nq = 3.6 x 10^"^ and To = 1.5 x 10^ yrs, and the disk's 
equilibrium timescale is Tgg = 2.0 x 10^ yrs. 
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Fig. 14. — Normal optical depth r„(r) that results in the surface brightness profile shown in 
Figure [T3| assuming Qg = 0.1 and Qs = 1. The dust production rate for the Qs = I model 
is 300 times smaller than the = 0.1 model. 

It should be noted that the planetesimal disk's inferred erosion rate is very sensitive 
to the dust grain's optical properties. Equation fl3T]) shows that an edge-on disk's surface 
brightness is varies as i? oc BqNqR^^^ where Bq oc Qs, Nq oc p^Qrp^'^ , and Rmin oc Qrp/ P 
(Equations [3| [6| l28a[ and [30l) . which means that the inferred dust production rate varies as 
Md oc /q where fg = Qrp/Ql and Qrp = 1 — gQs- (Interestingly, the inferred dust production 
rate is insensitive to the assumed grain density p, which cancels out.) So if the dust grains are 
instead dark and have a scattering efficiency Qs = 0.1, then the radiation pressure efficiency 
is Qrp — 1 and fg ~ 100, which means that the dust production rate must be ~ 100 faster 
if the dust grains are dark. In this case, (3 Pic would have lost ~ Mdti, ~ 10^ earth-masses 
over its lifetime, which is implausible. Rather, it is more likely that the dust grains at /3 Pic 
are bright. 

The optical depth profile r„(r) that is inferred from (3 Pic's observed surface brightness 
B{x) is shown in Figure [TH Note that B{x) oc J QsTndi where the integration is along the 
line of sight, so the disk's optical depth is uncertain by a factor of Qs- The Qs = I curve 
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obtained here is similar to the optical depth profile that lAhmic et al.l ( 120091 ) inferred from 
the same dataset§|. The range of possible values for the total dust mass and total dust 
cross section are 4.7 < < 170 lunar masses and 1.2 x 10^° < < 1.6 x 10^^ km^, 
with the lower values for a disk composed of bright Qs = 1 dust grains, and the higher values 
for dark Qs = 0.1 dust. 



3.1. sensitivity to parameters 

The following explores the sensitivity of the above results for /3 Pic to the assumed 
model parameters. 

First, the dust produced by the simulated planetesimal disk has a. q = 2.5 size distri- 
bution, so most of the just-produced dust grains' mass is in the largest grains that actually 
contribute little to the debris disk's optical depth. Consequently, the dust production rate 
inferred above is somewhat sensitive to /3mm; which is the size parameter of the largest grains 
that are allowed in the model. The equilibrium model shown in Figure [13] has Pmin = 0.065 
(which corresponds to a maximum grain radius Rmax = 23 /im), and it requires a dust pro- 
duction rate of Md = 1.7 x 10^^ gm/sec when the model is fit to the the disk's observed 
surface brightness profile assuming Qg = 0.7. However, cutting off the dust size distribu- 
tion at Prnin = 0.1 and then optimizing the other parameters then results in a fit that is 
only marginally acceptable; that model has Rmax = 14 /im and Ma = 1.2 x 10^^ gm/sec. 
Nonetheless, models having f3min ^ 0.2 are ruled out, because they are deficient in grains 
having sizes R > 6 fim that would be confined to radial distances of 75 < r < 250 AU, and 
their absence from the model reduces the simulated disk's surface brightness at projected 
distances a: ~ r. 

Interestingly, the disk's surface brightness profile B{x) is rather sensitive to the size 
distribution power law q. For instance, models having q < 2.0 do not reproduce the disk's 
observed surface density profile, because the model disk's inner r < 150 AU region is over- 
dense with large dust grains that contribute too much surface brightness at a; < 150 AU, 
while models having q > 3.0 are underdense in large grains and thus too dim there. Evi- 
dently, the dust being produced by the unseen planetesimals at f3 Pic have a g ~ 2.5 size 
distribution that is just a bit shallower than the canonical Dohnanyi q = 3.5 size distribution. 
And because the observed dust size distribution is constrained to 2 < g < 3, this translates 
into an uncertainty in the value of quoted above is no more than a factor of 2. Also, do 
keep in mind that the value for q refers to the size distribution for dust production by the 



Ahmic et al. 



(120091) work does not mention Q^; we suspect that Qs — 1 was assumed there. 
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planetesimal disk, and that the debris disk's resulting dust size distribution is actually much 
steeper due to collisions among dust grains (e.g., Figure [3]). 

This kind of modeling is also sensitive to the planetesimal disk's inner and outer radii. 
The midpoint of the knee at x ^ 110 AU in /3 Pic's surface brightness profile indicates that 
the planetesimal disk's midpoint is near r ~ 110 AU (Figure [T3l) . and best agreement with 
the observations is achieved when — 75 AU and Vout — 150 AU i.e., when the disk has a 
radial aspect ratio of Tin/rout — 0.5. However, broader disks having Tin/ Tout < 0.3 are ruled 
out because the knee in their surface brightness profiles are too broad to fit the observations. 
Likewise, narrower planetesimal disks having Tin/ Tout > 0.8 are ruled out because those 
models produce a prominent bump in the debris disk's surface brightness where x = r^^, and 
such bumps are not seen in these optical observa tions. Although these surface brightnes s 



bumps are seen in infrared observations of /3 Pic (IWahhaj et al.ll2003l : iTelesco et al.ll2005l ). 
their locations vary with the observing wavelength, which indicates that their interpretation 
likely requires modeling the disk's thermal emission as well. 

We have also considered 'inside-out' erosion of the planetesimal disk, with c = —2, but 
the edge-on disk's surface brightness was rather similar to that generated by planetesimal 
disk suffering 'outside-in' erosion with c = +2, so c = is adopted in all models shown here. 
Also note that strong dust grains having Q* = 10® ergs/gm are preferred over models that 
use weaker Q* < 10^ ergs/gm dust grains. The larger grain strength increases the survival 
of the larger dust grains (see Figure [3]) that are confined to the vicinity of the planetesimal 
disk at 75 < r < 150 AU. When Q* < 10^ ergs/gm, the lower abundance of large dust grains 
then results in a noticeable surface brightness deficit at x ~ 100 AU. 



3.2. implications for planet formation 

If one adopts the plausible assumption that the (3 Pic dust grains are bright {Qs ~ 0.7), 
then the mass of the planetesimal disk is likely well above 110 M^, which is sufficient to 
form about six Neptune-class planets. This planetesimal disk's considerable mass, plus its 
broad radial extent (75 < t < 150 AU), also provides several interesting constraints on the 
prospects for planet formation at (3 Pic. Analogy with the Kuiper Belt, which resides at 
distances ~ 30% beyond Neptune's orbit, suggests that there should be no Neptune-class 
giant planets orbiting well beyond r > 60 AU from /3 Pic, because such planets would have 
scattered away a portion of the inferred planetesimal disk and part of its dust disk. The fact 
that there are no such planets yet means that if planets are trying to form there via core- 
accretion (e.g., a runaway accretion of planetesimals by a growing planetary embryo), then 
the timescale for this kind of assembly in the 75 < r < 150 AU zone must be longer than the 
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system's ~ 12 Myr age. However, gravitational instability (GI) is a much faster planet- 
formation mechanism, one that can also form gian t planets at large stellarcentric distances 



(jPurisen et al.l 120071 : iDodson- Robinson et al.ll2009l ). However, GIs are only operative when 



the circumstellar disk is still rich in nebula gas, and the lifetime of a gas disk is at most a 
few Myrs. Consequently, the presence of (3 Pic's broad and massive planetesimal disk also 
tells us that planet formation via GI did not happen in the 75 < r < 150 AU zone. 

It can also be concluded that giant planets are unlikely to have formed interior to r ~ 75 
AU and subsequently migrated or scattered deep into the 75 < r < 150 AU zone. For if that 
were to have happened at /3 Pic, then the migrating/scattering planet would have accreted 
and/or scattered the planetesimal disk, causing the circumstellar dust to transition into a 
relic debris disk within a time ttrans ~ O.lTeg ~ 2 Myrs. Consequently, the probability of 
an astronomer observing this system before its surface brightness transitione d into a relic 
B(x) PC profile is ttrans/U ~ 15%. This finding is also consistent with iBooth et al. 



( I2OO9I ). who show that zero of 106 stars having circumstellar dust also appear to have expe- 



rienced a Nice-model scattering of giant planets. 

Of course, the Nice model was not intended to explain observations of circumstellar 
debris disks. However, our purpose here is to show how disk observations can be used to 
determine whether any particular Solar System's scenario might also be ongoing at other 
circumstellar disks, or if that scenario should be regarded as a rarity among the known 
population of star-disk systems. Indeed, it is a curiosity that (3 Pic seems to be planetless at 
r > 75 AU zone, despite having a substantial amount of planetesimal mass there. However 
(3 Pic's heavy dust production rate, which is at least ~ 9 M®/Myr (and may be much 
larger if the grains are darker than Qs = 0.7), suggests an alternate interpretation — that the 
planetesimal disk orbiting (3 Pic is not a planet-producing disk, but is instead a planetesimal- 
destroying disk, due to a vigorous collisional grinding of planetesimals and a heavy mass-loss 
rate due to the blowout of dust by radiation pressure. 



3.3. Testing the assumptions 

This subsection examines the key assumptions that are employed here. Those assump- 
tions are: (1) that collision fragments contribute little to the disk's optical depth, (2) the 
dusty disk is so tenuous that the dust grains are not shadowing each other, (3) that Poynt- 
ing Robertson drag is insignificant, and (4) that the planetesimal disk's dust production is 
steady over time. A final subsection then comments on the fact that many debris disks are 
non-axisymmetric. 



3.3.1. dust fragmentation 



When dust grains collide, they spawn smaller dust fragments, but their contribution 
to the debris disk's optical depth is ignored by this model. This is appropriate if the bulk 
of those fragments are sufficiently small and fast such that they are unbound and leave 
the system. It turns out that this assumption is quite reliable, and is confirmed from the 
model output. Begin with a single term in Equation (|9]), IZijk = a^kninj/Tout, which is the 
rate at which the dust in the target streamline i are destroyed due to collisions with the 
dust in streamline j at the site r = r^-fc where their orbits intersect. The resulting dust 
fragments are presumed to have a power-law size distribution Nf = CRJ'^, which gives the 
number of dust fragments of radius Rf produced by the destruction of a single target grain 
of radius Ri, with the coefficient C determined by mass conservation. Equations f lT^] - [T5|) 
then provide the speed vj at which the fragments recoil from the collision, and their specific 
energy E = ^vj — GM^{1 — f3f)/rijk indicates whether the collision fragments having a size 
parameter Pf will stay bound to the star or leave the system. The total rate at which 
the debris disk produces collision fragments of radius Rf is TZf = NfTZ^k where the sum 
proceeds over all sites in the disk rijk where dust grains collide and produce bound fragments 
of radii Rj. This is the secondary dust production rate, and it is to be compared to the 
planetesimal disk's dust production rate Pf. The ratio TZf/Pf depends of the radius of the 
smallest possible fragment Rg. Also recall that Rmin is the radius of the smallest bound dust 
grain that is produced in the planetesimal disk, so Rs < Rmin- In this instance, the above 
algorithm finds that TZf ^ Pf, which means that the rate at which collision are injecting 
second-generation dust fragments into the debris disk is negligible in comparison to the rate 
at which the planetesimal disk is injecting first-generation dust into the disk. For example, 
when Rs < O.dRmin, this algorithm finds that TZj/Pf < 1% when applied to the example 
model of Section [2.4.31 The ratio TZf/Pf is small because most of the collision fragments are 
small, unbound, and do not contribute to the debris disk, which justifies the model's neglect 
of dust fragmentation. 



3.3.2. shadowing 



If the disk is sufficiently dense, then dust grains can shadow each other and alter the 
disk's surface brightness profile B{x). The calculations presented above all assumed that the 
dust grains are not shadowing eac h other, and that assu mption is confirmed by examining 
the disk's radial optical depth ( jWeidenschillingI l2010[ l. Begin with a narrow annulus in 
the disk of radius r and radial width Ar. The cross section of dust in that annulus is A A = 
Tn{r)2'Kr Ar where Tn{r) is the disk's normal optical depth. The annulus' vertical thickness 
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is 2/r, so its optical depth along a radial line-of-sight is Axr = AA/27rr2/r = r„Ar/2/r. 
The total radial optical depth of dust interior to r is then 

..M^r^. (33) 

Jo 2Ir' 

which is easily integrated numerically. When this quantity is not small, then the dust grains 
at r are shadowed by interior dust, which reduces the stellar illumination there by a factor 
e"^''. But if ^ 1 then the dust grains are fully illuminated. 

Most of the debris disk models described in the previous sections have < 10~^. The 
two exceptions are the heavy dust-producing models, such as those having Md = 10^^ gm/sec 
that have ~ 0.01 (Figures HI O and|9]), while the (3 Pic model that considers dark Qs = 0.1 
dust grains has ~ 0.04. So in summary, shadowing by dust grains is not significant for 
the debris-disk models considered here. 



3.3.3. Poynting Robertson drag 

Poynting- Robertson (PR) drag is a weak force that causes the orbits of small dust grains 
to slowly decay. The following shows that this orbital decay is neghgible for the dust in the 
(3 Pic debris disk. 

PR drag is the acceleration that results when a moving grain absorbs and/or scatters 
stellar photons. Because a grain's motion also results in a slight 'headwind' of photons, 
the transfer of momentum, from the photons t o the grain, cause s its orbit to decay. The 



acceleration on a dust grain due to PR drag is ( iBurns et al.l 119791 ) 



apfl = ^ — r H e (34) 

\ c c J 

where Vr and vg are the grain's radial and tangential velocities (Equations [T3|) and f , 6 
are unit vectors in a polar coordinate system. Inserting this into the Lagrange planetary 
equations and time- averaging those equations over an orbit then provides the rates at which 
an orbiting dust grain's semimajor axis a and eccentricity e decay due to PR drag. 



ac 



[1 - e2)3/2 



and e = * — (35b) 



( IWvatt fc Whipple! 



19501 : 



Burns et al 



19791 ). The timescales associated with this orbit decay 
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are 



Ta = 1^ = ' TpR 36a 



e A^l^ I a 



and Te = —7 = ^ TpR (36b) 

|e| 5/3 V^o«t/ 



2 



where the constant 



Tp^ = ik^ = 2.0 X 10^ ( ^] ( -^^) ' yrs. (37) 
2GM^ \MqJ V50AU/ ^ ^ ' 

Inserting Equations ([5]) and ([7]) into the above shows that these timescales are simple func- 
tions of grain size /3 or R. If these orbit decay timescales are long compared to a dust grains' 
collisional lifetime Tc, then it is appropriate to ignore PR drag. 

These orbit decay timescale are plotted in Figure H] for the example debris disk that 
is generated by a planetesimal ring of radius rout = 50 AU (see Sections 12.4.31 - 12.4.51 and 
Figures [21 4TT]) . All of these curves diverge for small grains that have radii R ~ Rmin and size 
parameters /3 ~ |. Note that for most grain sizes Tc{R) <^ Ta{R), which means that the 
effects of PR drag are negligible over the lifetime of grains that are only slightly larger than 
Rmin- Figure m also shows shows that PR drag is only significant when the system's dust 
production rate Md is sufficiently low, which can make > large, but only for these small 
R ~ Rmin grains. But this only occurs when Md S. 10 ^ ^ gm/ sec, for the scenario considered 
in Figure m A similar result was obtained by IWyattI (120051 ) . who showed that PR drag is 
only significant when the debris disk is sufficiently tenuous. And finally, note that PR drag 
is completely negligible (T^ <^ T^) for grains of all sizes in the (3 Pic debris disk, due to its 
very vigorous dust production rate and the grains' short collisional lifetimes. 



3.3.4- dust production over time 

The model employed here assumes that the planetesimal disk's dust production rate 
is steady over time. But this assumption might seem debatable, because the collisional 
erosion that drives dust production will ultimately decrease the planetesimal disk's mass 
over time. But keep in mind that accretion within the planetesimal disk tends to produce 
larger bodies whose gravity can stir-up the planetesimal disk, and that can instead increase 
the disk's dust production rate. Further, simulations of the coUisional/accretional evolution 
of a planetesimal disk shows that the outc omes are very sensitiv e to the planetesimal disk's 



initial conditions that are poorly known (jWeidenschillingI |2010[ ). Due to this uncertainty 



and that any rigorous treatment of the planetesimal disk's evolution goes well beyond the 
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intended scope of this study, this model makes the simplest possible assumption, that the 
planetesimal disk's dust production r ate is steady over time. This assump t ion is also justified 
by the collisional evolution models of IStern fc Colwelll (Il997l ) and iKenyonI (j2002[ ). who obtain 
a slow erosion timescale of ~ 500 Myrs for a Kuiper Belt orbiting at r ~ 40 AU. Because the 
timescale over which the erosion-rate varies is likely longer than f3 Pic's equilibrium timescale 
(Teg ~ 20 Myrs, Section [3]), the debris disk is expected to remain in quasi-static equilibrium 
as it adjusts to any slow change in the dust production rate. Consequently, the debris-disk 
mass and dust production rate inferred here are expected to be reliable. 

Also note that /3 Pic's dust production rate is considerable, at least Md ~ QMg/Myr. 
Consequently, steady dust production over the age of the system also requires a large 
reservoir of planetesimal mass Mp, with Mp ^ M^^tj, ~ 110 M^, i.e., the mass in /3 Pic's 
Kuiper Belt at 75 < r < 150 AU must be well in excess of 110 earth masses. Alternatively, 
if /3 Pic's Kuiper Belt does not satisfy Mp ^110 M^, then collisions among planetesimals 
plus blowout by radiation pressure is going to grind down this extra-solar Kuiper Belt in a 
fewxtj, ~ lO's of Myr. In this case the assumption of steady dust production does not hold, 
but it does mean that dust production was probably more vigorous in the past. Regardless, 
P Pic's planetesimal disk is or was very massive. 



3.3.5. disk asymmetry 



Many edge-on debris disk s are lopsided, wi th one ans a substantially b righter than the 
other; examples in clude f3 Pic JHeap et allboooh . AU Mic JKrist et al.ll2005[ l. and HD 15115 
(IKalas et al.ll2007l ). It may be that circumstellar debris disks are routinely non-axisymmetric. 
Once possible explanation for a disk's asymmetry is a recent dust-producing collision in the 
planetesimal disk. Recall that the orbits of all the dust grains produced in a single collision 
have their longitudes of periapse u aligned (Section [2TT]) . Because dust produced in a collision 
tend to loiter at their apoapse, the debris disk will have an excess of dust at longitudes 
6 = u + 7:. So if that collision was vigorous enough to produce a substantial amount of dust, 
then the disk will be non-axisymmetric, and may appear lopsided when viewed edge on. 



Another possible explanation for a disk's asymmetry is described in iHahnl (120091 ). who 
showed that if the circumstellar dust is produced by an eccentric planetesimal ring, then the 
smaller dust grains produced at the ring's periapse are less bound to the star due to their 
higher orbital velocities there. Consequently, there will be fewer small and marginally bound 
grains in the direction of the ring's apoapse, which will make that part of the debris disk 
dimmer than the periapse side that also has an excess of small grains. This is a particularly 
interesting scenario since it also implies that an extra-solar planet is likely present in order 
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to maintain the planetesimal ring's forced eccentricity. 

Both of these scenarios will be explored in a foUowup study using a more advanced 
version of the debris-disk model that will be generalized to handle the coUisional evolution 
of a non-axisymmetric debris disk. 

Lastly, it should be noted that this model also ignores the planetesimals' free eccentricity 
e f that is associat e d with the random part of their noncircular motions. Equations (13-14) 



m 



Thebault et al.l ( 120031 ) show that the planetesimals' random motions tends to blur the 
simple relationship between a dust grain's orbit and its size parameter (e.g., Equations [5] 
and[7j). However this effect will be small when the planetesimals eccentricity e/ is small, 
which in fact is a requirement in order for planetesimals to have formed in the first place. 



4. Summary and Conclusions 

A numerical model for a circumstellar debris disk is developed and applied to observa- 
tions of P Pictoris. The model accounts for dust production by colliding planetesimals and 
dust destruction due to collisions among grains. These rates for dust production and coUi- 
sional destruction also provide a rate equation whose solution gives the dust abundance over 
time and as a function of grain size. That solution also provides the debris disk's grain size 
distribution, which is steepened substantially by collisions among dust grains. These calcula- 
tions also give the dust grains' coUisional lifetime Tc{R), which depends on the planetesimal 
disk's dust production rate Md and grain radius R. 

Scaling laws are derived, and it is shown that the dust abundance N in the debris disk 

• ]^ /2 

varies as oc once the disk has settled into coUisional equilibrium, and that dust 

■ —1/2 

lifetimes vary as Tc oc . It is also shown that the radial drift of dust grains due to PR 

drag is unimportant provided the dust production rate is sufficiently high so that the grains' 
coUisional lifetimes are short comp ared to the drift timescale. The model also recovers the 



results of IStrubbe fc Chiangl (120061 ). who showed that the debris disk's optical depth varies 
as Tn{r) oc r~'^/^ when in equilibrium, and that an edge-on disk's surface brightness varies 
as B(x) oc where r and x are the radial and projected distances from the central 

star. Note though that those optical depth and surface brightness profiles will be steeper if 
the disk is young and not yet in coUisional equilibrium. Alternatively, those profiles will be 
shallower if dust production has ceased in the past, which might occur if planets form in (or 
migrate or scatter into) the planetesimal disk and cause its dynamical depletion. 

It is also shown that these quantities are not very sensitive to the dust grains' strength 
Q* when Q* < 10^ ergs/gm. However, an edge-on disk's surface brightness profile is quite 
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sensitive to the dust grain's light scattering asymmetry parameter g, and disks hke the ones 
at P Pic and AU Mic having a 'knee' in their surface brightness indicates that the dust are 
very asymmetric hght scatters, with \g\ > 0.7. 

The model's principal dynamical parameters are the planetesimal disk's radius and its 
dust production rate, so a comparison to observations then yields estimates of or else limits 
on these important cosmogonic quantities. For instance, fitting the model to optical HST 
observations of (3 Pic shows good agreement with the disk's observed surface brightness 
profile when the unseen dust-producing planetesimal disk there is quite broad, extending 
over 75 < r < 150 AU. This disk's dust production rate is also quite heavy, Md ~ 9 
M^/Myr, if it is assumed that the dust grains are bright like Saturn's rings {Qs = 0.7). 
In this case, the total cross section of dust in the debris disk is = 1.9 x 10^° km^, and 
the inferred mass there is =11 lunar masses. Note that this m ass is comparable t o 



that previously inferred from submillimeter observations of this disk (IHoUand et al.lll998l ). 
Also, it is unlikely that the dust grains at (3 Pic are dark {i.e., with Qs < 0.1), since that 
would require the planetesimal disk to grind away at a rate that is implausible. Indeed, the 
mass-loss rates inferred here are so heavy as to suggest that the r > 75 AU zone at P Pic 
might be a region of planetesimal destruction due to collisional grinding, rather than a site 
of ongoing planet formation. 

The model developed here is called ddisk, which is an easy-to-use IDL script that others 
might wish to use as they diagnose their observations of circumstellar debris disks. This code 
is available for download at http : //gemelli . spacescience . org/ ^haJinjm/ software . htmlj 
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